A=[5 1/3; 1/3 1];
syms x1 x2
x=[x1 x2].';
A*x
simplify(x.'*A*x)
x=[6 1]';
x'*A*x
x=[1 3]';
x'*A*x
